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We introduce a strongly correlated mean-field theory of the ionic Hubbard model on the triangular 
lattice with alternating stripes of site energy using Barnes- Coleman slave bosons. We study the 
paramagnetic phases of this theory at three quarters filling, where it is a model of Nao.sCoC^, 
Rbo.5Co02, and K0.5C0O2. This theory has two bands of fermionic quasi-particles: one of which 
is filled or nearly filled and hence weakly correlated; the other is half- filled or nearly half-filled and 
hence strongly correlated. Further results depend strongly on the sign of the hopping integral, t. 
The light band is always filled for t > 0, but only becomes filled for \A/t\ > 1.5 for t < 0, where 
A is the difference in the site energies of the two sublattices. A metal-charge transfer insulator 
transition occurs at \A/t\ = 5.0 for t > and |A/t| = 8.0 for t < 0. In the charge transfer insulator 
complete charge disproportionation occurs: one sublattice is filled and the other in half filled. 
We compare our results with exact diagonalisation calculations and experiments on Nao.sCoC^, 
and discuss the relevance of our results to Rbo.sCoC>2 and K0.5C0O2. We propose a resolution of 
seemingly contradictory experimental results on Nao.sCoCb. Many experiments suggest that there 
is a charge gap, yet quantum oscillations are observed suggesting the existence of quasiparticle states 
at arbitrarily low excitation energies. We argue that the heavy band is gapped while the light band, 
which contains less than 1 charge carrier per 100 unit cells, remains ungapped. 

PACS numbers: 



I. INTRODUCTION 

Experiments on A^CoC^, where A is Na, Rb, or 
K, show a wide range of strongly correlated phases* 1 * 2 - 
Nao.5Co02, in particular, has attracted much attention 
following the discovery of superconductivity when water 
is intercalated into the system^ The phase diagram, with 
doping x, shows many interesting phases including 
a 'Curie- Weiss metal', A type antiferromagnetism (fer- 
romagnetic layers stacked antiferromagnetically) and an 
'insulating' phase seen only at x = 0.5. The latter 
phase is particularly puzzling as many probes [includ- 
ing resistivity^ optical conductivity^ and angle resolved 
photoemission spectroscopy (ARPES)^] suggest that it 
is insulating, however, Shubnikov-de Haas oscillations are 
also observed^ suggesting the state is metallic. The aims 
of this paper are to present a simple variational theory 
of j4o.5CoC>2, where A=Na, Rb, or K, and to attempt to 
reconcile these seemingly contradictory experiments. 

In Aj;Co02 the Co atoms form a triangular lattice 
and the simplest model of the band structure, a single 
band triangular lattice with nearest neighbour hopping 
only, gives good agreement with ARPES experiments on 
Na^CoC^^ The doping of the system is controlled by 
the concentration of A ions, x, with the single band be- 
ing half filled at x = and filled at x = 1. However, 
the A ions order so as to minimise the mutual Coulomb 
rcpulsion i 9 ' 10 : 11 Therefore, the Coulomb potential due to 
the A ions is different at different Co atoms and, because 
the A ions are ordered, this gives rise to an ordered ar- 
rangement of potentials at the vertices of the triangular 



lattice. Thus an effective Hamiltonian for ^CoO^ is the 
ionic Hubbard models 

In addition to the interest in A^CoU^ the ionic Hub- 
bard model is of significant interest in its own right. The 
ionic Hubbard model on half filled bipartite lattices has 
attracted interest because it undergoes a transition from 
a Mott insulator to a band irismatori2iIiyiiI&iI2iILI£ 
Furthermore, away from half filling and on frustrated lat- 
tices the ionic Hubbard model shows a subtle interplay 
between charge and spin ordering and metallic and insu- 
lating phases i 19 i 20 ' 21 i 22 ' 23 i 24 The Hamiltonian of the ionic 
Hubbard model is 

(ij)a- i ia 

where t is the hopping amplitude between nearest neigh- 
bour sites only, U is the effective on-site Coulomb re- 
pulsion between two electrons, e, is the site energy, cfj 
annihilates (creates) an electron on site i with spin a 

and ilia — cj^Cjo-. Previous studies of this model on 
frustrated lattices included both analytical and numer- 
ical studies of zigzag ladders ; 19 ' 20 dynamical mean field 
theory (DMFT) studies on infinite dimensional fco^i and 
two-dimensional triangular lattices^ and exact diagonal- 
isation on small triangular lattices! 23 ' 24 A more extensive 
discussion of previous work on the ionic Hubbard model 
is given in Ref. 24;. However, here we aim to provide a 
simple variational description that captures as much of 
the strongly correlated physics of this model as possible. 

The wide range of numerical techniques, described 
above, that have been applied to the ionic Hubbard 



model have not, previously, been complemented by a 
commensurate effort to develop simple variational ap- 
proaches. The history of the theory of strongly correlated 
electrons shows that progress has often been made when 
accurate numerical techniques are combined with such 
variational calculations. Therefore, our theory, which 
can be straightforwardly generalised to other potential 
arrangements, lattices, fillings, etc., is complementary to 
the previous numerical work. 

In Na^OoC^ many different Na ordering patterns are 
seen at different values of cc^ i 10 ? 11 and each of these corre- 
spond to a different {fii}^ In principle, each of these dif- 
ferent {ei} correspond to a different Hamiltonian. There- 
fore, for simplicity and definiteness, we specialise to the 
Hamiltonian relevant to x — 0.5. This model has two 
sublattices with different site energies, e, = A/ 2 for A- 
sites and t\ = —A/2 for B-sites,— and with the sublattices 
arranged in stripes as shown in Fig. [TJs. Thus the Hamil- 
tonian is 

U = -t (cjaaCj/Sff) +Uj2 fl ^n tai 

{ija/3)a ia 

+ ~2 ( r " lL4<J ~~ hiBa) , (2) 
ia 

where c^ CT annihilates (creates) an electron with spin a in 
an orbital centred on site i belonging to sublattice a and 
(ija(3) indicates that the sum is over nearest neighbours 
on the appropriate sublattices only. 

Hamiltonian (|2j) is a natural instance of Hamiltonian 
([T]) for initial study on both theoretical and experimental 
grounds. Experimentally, the behaviour of, particularly, 
Nao.5Co02 is very different from that of Na^OoC^ with 
x either a little larger or smaller than 0.5. Further, 
many of the phenomena observed at x = 0.5 seem to 
have a natural explanation in terms of the ionic Hubbard 
model.— i 23 ' 24 We have previously studied the quarter- 
filled ionic Hubbard model with a stripe potential by ex- 
actly diagonalising small clusters. 23 ! 24 We found a com- 
plicated interplay between the charge and spin degrees 
of freedom and between insulating and metallic states. 
We also found that these calculations provide a possible 
framework for understanding a wide variety of experi- 
ments on Nao.5Co02, Rbo.sCoC^, and K0.5C0O2. How- 
ever, the most interesting regime (i.e., — t ~ A -C U, 
see section ITBj) . in terms of its relevance to experiments, 
is also the most challenging to investigate theoretically. 
This provides additional motivation to investigate the 
simpler mean-field theory presented below. Comparison 
between our mean-field theory and these exact diagonal- 
isation calculations will be made, where possible, below. 

The rest of this paper is organised as follows. In the 
remainder of this introduction we briefly review the most 
pertinent experiments on Nao.sCoC^ and previous theo- 
ries aimed at explaining these measurements. In section 
im we give the formal details of our theory and derive a 
set of self-consistency equations for the ground state. In 
section IIIII we report the results of numerical solutions 




FIG. 1: [Color online] The various unit cells and Brillouin 
zones discussed in this paper: (a) the unit cell of the isotropic 
triangular lattice; (b) the unit cell of the ionic Hubbard model 
[Eq. @]; and (c) the basal plane of the unit cell of A0.5C0O2; 
panels (d), (e), and (f) show the Brillouin zones corresponding 
to the unit cells in panels (a), (b), and (c) respectively. In (b) 
the different size circles distinguish the A and B sublattices. 
In (c) the large red (small blue) circles indicate the positions 
of Na ions above (below) the C0O2 plane. The unit cell of 
the ionic Hubbard model (b) is twice as large as that of the 
Hubbard model of the triangular lattice (a) . Further note that 
the basal plane of the unit cell of ^4o.sCo02 is twice as large as 
the unit cell of the ionic Hubbard because of the periodicity 
of the (out of plane) A ions. Selected high symmetry points 
are marked in the Brillouin zones. Note that the Brillouin 
zones are drawn to scale and that the F-point [k = (0, 0)] is 
equivalent in each Brillouin zone. 



of these equations and compare our results with exper- 
iments on 74.0.5C0O2 and previous theoretical studies of 
the triangular lattice ionic Hubbard model. In section 
IIVI we draw our conclusions. In appendix lAl we describe 
the band structure of the model © for non-interacting 
electrons. 



A. Experiments on Nao.r,Co02 

There has been far less experimental work on 
K0.5C0O2 or Rbo.5Co02 than Nao.sCoC^. Therefore, 
in this section, we focus mainly on experiments on 
Nao.5Co02- The brief review below brings out two par- 
ticular puzzles that theory needs to explain about the 
ground state. The first is how it is some properties are 
consistent with a metallic ground state and that oth- 
ers are more consistent with an insulating ground state. 
The second puzzles concerns how it is that the ground 
state has a large magnetic moment but a small amount 
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of charge ordering. 

1. Charge order on the cobalt ions 

Measurements of the relaxation rates for 59 Co NMR 
found two distinct Co sites, consistent with a charge or- 
dered stated However, the authors noted that the degree 
of charge disproportionation is rather small. More recent 
59 Co and 23 Na NMR measurements find no signature of 
differences in the charge state of the cobalt ions at the 
two distinct sites^ suggesting that all the Co atoms are 
in the charge state Co( 3 ' 5±1 " + , where 77 < 0.2 gives an 
upper bound on the extent of charge ordering^ This is 
in contrast, to the case of x = 0.7 or 0.75 for which NMR 
clearly detecting a charge disproportionationf 3 ^- 35 High 
resolution neutron crystallography detects small differ- 
ences in Co-0 bond lengths, and an analysis based on 
bond valence sums, is consistent with a charge order of 
77 ~ 0.06. 27 Very recent NMR measurements were able 
to detect a very small charge ordering along the Co(l) 
chains, below 51 Kj 2 ^ but the authors do not discuss the 
magnitude of the difference in charge between two sub- 
lattices. 



2. Transport properties 

As the temperature decreases the intralayer resistivity 
increases monotonically. No feature is seen at the mag- 
netic ordering temperature (88 K, see below). Above 
51 K, it is weakly temperature dependent with values 
of a few mfi cm,— characteristic of a "bad" or incoher- 
ent metal which does not have well-defined quasi-particle 
excitations* 2 - Below 51 K, the resistivity increases signifi- 
cantly, consistent with an activated form with energy gap 
of about 10 meV4 Although, it is usually stated that the 
transition at 51 K is a "metal-insulator" transition, we 
stress that it is really above 51 K Nao.sCoO^ is a bad 
metal and below 51 K the experimental evidence is not 
all consistent with the claim that Nao.sCoO^ is an insu- 
lator (see section ITA4[) . 

Applying hydrostatic pressure causes the resistivity to 
decrease and the temperature where the resistivity be- 
comes activated to increase. Above about 13 GPa the 
resistivity has a metallic temperature dependence^ In 
contrast, for x—0.75 Naa;Co02 is metallic at ambient 
pressure but becomes insulating above about 23 GPa. 

Resistivity measurements found that the charge gap 
apparent from the resistivity was suppressed by magnetic 
fields parallel to the layers, larger than about 35 T, with 
hysteresis between 15 and 40 T, at low temperatures* 7 - 
A field perpendicular to the layers does not suppress the 
"insulating" stated 

The thermopower and Hall coefficient have small and 
weakly temperature dependent positive values above 100 
However, they change sign near 88 K, obtain large 
negative values below 50 K, and then decrease towards 



zero at low temperatures. The latter is distinctly differ- 
ent from the behaviour of conventional semiconductors 
and insulators, for which the thermopower and Hall con- 
stant both diverge as the carrier density vanishes with 
decreasing temperature. 

3. Other evidence for a gap 

Quasi-particle features are only observed in angle 
resolved photoemission spectroscopy (ARPES) below 
about 150 K and reveal a Fermi surface consistent with a 
hole-like band for a tight-binding model on the triangular 
lattice with an effective hopping integral of t e tf ~ — 16 
meVpi about six times smaller than that predicted by 
LDA calculations. A gap begins to open up over the 
Fermi surface below about 60 K, and has a magnitude of 
about 6-11 meV at 10 K£ 

Measurements of the frequency dependence of the con- 
ductivity show no evidence of a Drude peak, consistent 
with the bad metal behaviour, seen in the static transport 
quantities.— £ Below about 100 K there is a loss of spectral 
weight at frequencies below about 100 cm -1 , consistent 
with the opening of a charge gap. This leads to a peak 
in cr(ui) at u ~ 200 cm -1 .- In comparison for Naa;Co02, 
with x away from 0.5, a Drude peak is present at low 
temperatures and the total spectral weight of the optical 
conductivity scales with 1 — x& 

4- Evidence for a metal 

Shubnikov-de Haas oscillations are observed in the field 
range of 15-30 Tf 7 - A fast Fourier transform of the oscilla- 
tory part of the conductivity has peaks at frequencies of 
150 T and 40 T. These frequencies correspond to pock- 
ets with cross-sectional areas of 0.25 % and 0.06 % of 
the area of the undistorted hexagonal Brillouin zone (cf. 
Figure Id). The effective mass associated with these or- 
bits are (in units of the free electron mass) 1.2 ± 0.1 and 
0.6±0.5, respectively. These measurements are extremely 
surprising given the range of other evidence suggesting 
that Nao.5Co02 has an insulating ground state. 

5. Magnetic properties 

Muon spin rotation (/iSR) experiments first saw evi- 
dence for magnetic order below about 86 K, and slight 
effects at 50 K and 30K. 33 Below 88 K, there is a splitting 
of NMR lines, consistent with the development of com- 
mensurate antiferromagnetic order The change in the 
resistivity from a bad metal to an activated behaviour 
at 51 K has little effect on the magnetic state. Elas- 
tic neutron scattering also detects long-range antiferro- 
magnetic order below 88 K, with a magnetic moment of 
0.26(2)/xb per magnetic cobalt ionj 31 i 37 The authors in- 
terpret this order in terms of alternating rows of magnetic 
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Co 3+ and non-magnetic Co 4+ ions, (i.e., complete charge 
order). However, we stress that the observation of such 
magnetic order does not require charge ordering to be 
present. Simple classical arguments suggest the magnetic 
moment should be less than (j,b(it-b — t^a)/^ where n v is 
the average number of electrons on the v — A, B sublat- 
tice. Hence, the charge disproportionation observed in 
the crystallography experiments discussed above, would 
imply a moment of about 0.06/zs, about one quarter of 
the observed value. 

Both the bulk magnetic susceptibility and the Knight 
shift are weakly temperature dependent, with a peak 
around 300 K, and a magnitude of about 5 xlO -3 
emu/molj 4 * 2 ^ This suggests an antiferromagnetic ex- 
change between magnetic ions of order a few hundred 
Kelvin. 



6. Related materials 

At other values of x in Na^CoC^, the ordering of 
sodium ions has been found to be important for the ex- 
istence of various ordered phases . 34 ' 35 ' 36 ' 38 ' 39 No insulat- 
ing state is seen in the corresponding misfit cobaltates, 40 
which supports the claim that Na-ordering is necessary 
for the insulating state. Furthermore, cooling a material 
at different rates found that the presence of sodium or- 
dering can drive an additional magnetic phase transition 
at x = 0.8 and 0.8541 



B. Theories of Na 5C0O2 

Choy et alM- considered an extended Hubbard model 
that included the Coulomb interaction between electrons 
on neighbouring sites, V, but neglected the Na-ordering 
and the 'ionic' term in Eq. $2§. They argued that V 
stabilises a charge ordered state and considered an ef- 
fective low-energy Hamiltonian for this charge ordered 
state. By fine tuning the parameters in their effective 
Hamiltonian Choy et al. were able to reproduce the ob- 
served temperature dependence of the resistivity and the 
Hall coefficient. 

Lee et a/ . 43 ' 44 have studied Nao.sCo02 via LDA+U cal- 
culations that allow for the effects of Na-ordering. They 
find a first order metal-charge ordered insulator transi- 
tion as U is increased. The charge disproportionation 
and the opening of the gap both occur at a single phase 
transition in this weak coupling theory. 

We have reported exact diagonalisation calculations 
for Hamiltonian @ of finite lattices. 23 ! 24 These studies 
suggested that, for parameters relevant to Nao.sCo02, 
the system is a covalent insulator, with little charge 
disproportionation and a small gap ~ |t|, rather than 
a charge ordered insulator, with strong charge dis- 
proportionation and a gap ~ A. This prediction is 
consistent with the experimental measurements of the 



charge disproportionatiom 26 ' 2741 Further, these calcula- 
tions predict a gap of the same order of magnitude as in 
seen in the ARPES 6 and the resistivity, 4 predict the large 
moment seen in neutron scattering ) 31 ' 37 and reproduce 
the main features seen in the optical conductivity^ Zhou 
and Wang 45 also studied a model that incorporates the 
effects of Na-ordering within the Gutzwiller approxima- 
tion and found it could explain many of the experimental 
observations. 

There has been significantly less work on 
K0.5C0O2 and RD0.5C0O2. Lee and Pickett 46 have 
reported LDA calculations for K0.5C0O2 and compared 
these with equivalent calculations for Nao.5Co02- 
They found that the t2 g band is rather narrower in 
K0.5C0O2 than Nao.5Co02 and that the Fermi surface 
is more complicated in K0.5C0O2 as there are several 
small pockets along the X-S line. They speculated that 
these pockets, some of which are quite well nested near 
the X point, may lead to enhanced magnetic tendencies 
in K0.5C0O2. 

Clearly an important question is: what parameters 
in Hamiltonian ((2|) correspond to the various possible 
choices of A in A0.5C0O2? Various atomistic calcula- 
tions that address this question have been presented for 
A=Na, but for other A there are not yet suitable esti- 
mates. CAS+DDCI calculations^ on small clusters of 
Nao. 5 Co0 2 give -t = 0.08 - 0.14 eV, U = 2.5 - 2.8 
eV, and A = 0.16 eV. For bulk Na . 5 CoO 2 the LDA 
yields A = 0.07 eV-2 and t ~ -0.1 eVi 3 -^, and elec- 
trostatic Ewald calculations give 1 ^ A = 0.03 — 0.05 eV. 
Furthermore, C0O2, which is described by Hamiltonian 
@ at half filling with A = 0, is observed to be a 
strongly correlated metal rather than a Mott insulator. 49 
This suggests that U < U c ~ 12—154, where U c is 
the critical value for formation of a Mott insulator on 
the triangular latticei 2 -^ Hence, realistic parameters for 
Nao.5Co02 may be in the range 10 < U/\t\ < 15, and 
0.5 < \A/t\ < 2. 

The negative sign of t is natural if the hopping between 
Co sites is dominated by the contribution from hopping 
via an intermediate O site. To leading order this gives 
t = -|4coo| 2 /(eco~eo), where t Co o is the direct hopping 
between a Co atom and a neighbouring O atom, eco is the 
energy of an electron in a d orbital centred on a Co atom, 
and eo is the energy of an electron in a p orbital centred 
on an O atom (note that one expects that eco > eo)- 



II. SLAVE BOSON THEORY OF 
STRONGLY-INTERACTING ELECTRONS 

We now investigate paramagnetic phases of the U — > 
00 limit of the ionic Hubbard model by deriving a 
mean field theory using Barnes-Coleman slave bosons^i 
The U — > 00 limit great simplifies the analysis and is 
not unreasonable in spite of our estimate, above, that 
10 < U/\t\ < 15. Our previous exact diagonalization 
calculation s 23 ' 24 show that, provided U \t\ the physics 
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of this model is insensitive to particular value of U/\t\. 
Specifically, our results for U — 10|t| show only minor 
quantitative differences from those for U = 100|t|. 

A. Theory of the metallic state 

We begin by making the particle-hole transformation 
Ciaa —> h\ au . We then introduce the slave bosons. Fol- 
lowing, for example, Ref. I52I we write 

hiacr = fiao-bia + e aa' fiacr' d\ a , (3) 
a' 

where e aa i is the completely antisymmetric tensor, 
annihilates (creates) a fermion with spin a at site i on 

sublattice a, d^J annihilates (creates) a boson, corre- 
sponding to a site that is doubly occupied by holes, at 
site i on sublattice a, and b^J annihilates (creates) a bo- 
son, corresponding to a site containing no holes, at site 



i on sublattice a. If the number of holes is less than or 
equal to the number of lattice sites and U = oo there is 
zero weight for configurations containing any sites dou- 
bly occupied (with holes) in the the ground state or any 
excited states at finite energy above the ground state. 
Therefore, d ia can be 'deleted' giving us 

hLv = fLbia- (4) 

It is clear from Eq. (TJJ that the f's are hole-like opera- 
tors. This transformation is exact, in the limit U — ► oo, 
provided we impose the constraint 

1 = Y / jl^ahfia + b ip bin = Qip (5) 
cr 

on the system. 

Hence, we find that the Hamiltonian ([2|) may be writ- 
ten as 



H-^N = t flJi^k* + y£ (b\ A k A - bisks) + (6) 

(ijaP)a i i/3 

We make a mean-field approximation by replacing the bosons by the expectation value q%j a p = (b\ a bj^), which 
gives, 

H mf - fj,N = t fLafiPrljiPa + y Y (quaa - QhBb) - H 53(1 + qupp) (7) 

(ijaf3)a i i/3 

We now assume that the bosonic mean field is homogeneous {qijap = Qap) an d introduce the Lagrange multipliers, 
A Q to enforce, on average, the local constraints [via the term — J~l ia A a (l — Qi a )]. Thus one finds that 

f = n w t - f iN-J2 Mi - 4a) (8) 

ict 

= 1 Y fLJjPoWa + y Y ( qAA - qbb) - vY^ 1 + iw) ~ Y x p ( 1 ~ q w ~ Y fipJip* ) • 

{ijaf3}cr i i/3 i/3 \ cr J 

Upon performing a Fourier transform one finds that there are two bands, which we denote the bonding (-) and 
antibonding (+) bands, 

P _ ( fkAcf \ f ( 2tq A A cos k x + X A Uq AB cos \ cos \ \ / f^Aa 
^ \ /kBcr / I 4:tq BA COS COS 2tq B B cos k x + \ B I \ fkBa 



+N l T ( qAA -Qbb)~Y M 1 + 9uu) + A„(l - q vu )] I 

Y 

kcr 



(9) 



V'k+cr V Z' Bk+CT ^ ( V'k+cr 

!>k-cr 



(10) 
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TABLE I: Physical meaning of the mean fields, qi, and the Lagrange multipliers, X±. Note that one can also interpret the 
quasiparticle weights as the inverse of the corresponding effective mass. 



Parameter Alternative symbol meaning 



q+ Z+ uniform intra-chain quasiparticle weight 

q- Z- anisotropic intra-chain quasiparticle weight 

|?ab| Zab inter-chain quasiparticle weight 

qAA Zaa A-sublattice intra-chain quasiparticle weight 

qBB Zbb B-sublattice intra-chain quasiparticle weight 

— A+ [i* effective quasihole chemical potential 

— A_ A*/2 effective ionic potential 



where k x and k y are defined in the reduced (1 x V3) Brillouin zone of the model (Fig. [T};), N = J^. is the number 
of unit cells, i.e., half the number of lattice sites, and the dispersion relations, E^± a , of the 'quasiholes', which are 
destroyed (created) by ^ffi.-, are given by 



/ 2 k x k 

i?k±cr = 2tq + cos k x +X + ±J [2fc?_ cos k x + A_] + 16t 2 \qAB\ 2 cos 2 cos 2 -y- 

where q + = ^ {qAA + qBB), q- = \ {qAA - qBB), X + = ^(A^ + X B ), and A_ = \{X A - A B ). 
Hence we find that 

F = ^ l 2t Q+ cosk x + A + + a\j [2tq- cos k x + A_] 2 + l6t 2 \q AB \ 2 cos 2 y cos 2 y 1 i>l a<7 i>k aa - 

ko<r V ) 



(11) 



-2N 



~+A_ ) 7- -/'(I + <1- ! - ■ A, CI -<m i 



(12) 



where a € {±1}- 

It is interesting to compare this with the solution of the model with U = [see Appendix [X] and particularly Eq. 
flAip ]. This enables a straightforward identification of the physical meaning of the mean fields q± and qAB and the 
Lagrange multipliers A± , see table HI 

It follows from the the approximation that the bosonic mean-field is homogeneous and the constraint of one particle 
per site that, at three-quarters filling, we need only consider a two site model (i.e., one unit cell), which must contain 
exactly one boson. Thus, the bosonic wavcfunction may be written as 



*>=I1 (ve»bl A + Vr^bt B ) |0> 



(13) 



where v is a real number between and 1. Thus 
qAA = (y\P iA bi A \V) = v 2 

qBB = (*|St B S iB |*) = l-v 2 

qAB = (tf |6t A S iB |tf) - v(l - v)e-* e = JqAAqBBer 16 
Therefore, the total energy of the model is 



F 



f) = e{ 2 ^ 

kacr k 



q\-q 2 _e^ 



cos k x + A4 



01/ [2tq^ cosk x + \_] + l&t 2 {q 2 + - q 2 _ ) cos 2 y cos 2 y \ n kao 



+2N 



+ A_ )q- - fl{l + q+) - A+(l - q+) 



(14a) 
(14b) 

(14c) 



(15) 



where n^ aty = (V>l acr ^kacr) • The phase 6 does not change the energy, thus, the solution has a U(l) degeneracy above 
that found in the slave boson mean field theory of the Hubbard models due to the two sublattice structure of the 
problem. 
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It is straightforward to show that q + = x = 1/2 follows from the requirement that —dF/d/j, = N e , where the total 
number of electrons is N e — J^i^a ^na^a = 2iV(l + x) for A2C0O2 and we have specialised to x = 1/2, which is 
the relevant filling. This result fits with our intuition that, on average, there should be one boson (i.e., one doubly 
occupied site) per unit cell. 

Similarly, (27V) -1 J2kaa n kaa = (1 — q+) = \ follows from the requirement that F is an extremum with respect to 
A+. Again, this fits with our intuition that there is, on average, one quasiholc fcrmion (singly occupied site) per unit 
cell. 

Three more complicated intertwined self-consistency conditions can also be derived: 

A 1 xr- 1 \2tq- cosk x + A J 2tcosk x - 16t 2 (?_ cos 2 %■ cos 2 % j 
1 ^ V kaa I J [2tq_ cos k x + A_] + 16t 2 (q 2 _ — g 2 ) cos 2 4f cos 2 



2 ^ UkJ 2 

follows from the requirement that F is a minimum with respect to q- ; 

1 ^ I 2tq- cos k x + A_ I 

1- = -^}^ a S 1 o ==>?w (16b) 

ZJV ka„ [ y [2tq_ cos fc x + A_] 2 + 16t 2 (q^ - q 2 _ ) cos 2 4f cos 2 J 

follows from the requirement that F is an extremum with respect to A_; and 

j 16f 2 g+ cos 2 cos 2 1 

A* = \ 2t COS Ka; + Ct > rikacr + A -1 



[2tq_ cos fc x + \-Y + m 2 {q\ - q 2 _) cos 2 \ cos 2 ^ 



follows from the requirement that F is a minimum with respect to q + . 

I 



(16c) 



In order to compare our results with ARPES exper- 
iments on Nao.5Co02 and K0.5C0O2 we calculate the 
spectral function, A(k, w) = - ImTrG fj, vlTlT < (k, cj), where 
G^vaa 1 (k, w) is the electronic propagator. In real space 
and time the one electron propagator is given by 

iG ijliVtTtT >(t) = (c^aWcjwCO)) 

= ( 17 ) 
Fourier transforming and taking the trace, one finds that 

A(k,u) = 2j2z ]ia [l-n F (E ]iaa )}6(E kaa -u;), (18) 

a 

where uf is the Fermi function, 

Z ka = q+ + a cos 28 k (19) 
is the quasiparticle weight, and 

tan fl k = £k+T = 2f (!ll^) COS fc - + A + (20) 
4tJq\ - q 2 _ cos %■ cos ^ 

Therefore, we find that, in contrast to moan field slave 
boson theories of other models, the quasiparticle weight is 
momentum dependent in the theory of the ionic Hubbard 
model. However, this does not mean that all non-local 
correlations induced by the on-site Coulomb repulsion 
are included in the present approach. 



B. Theory of the insulating state 

With any minimisation over a finite parameter space 
it is important to check explicitly whether any of the 
end points are the global minimal Therefore we must 
consider the special case g_ = ±1/2 separately from the 
general set of self consistency conditions [Eqs. fTo])]- Let 
us consider the case of large (how large will be determined 
later) positive A such that g_ = —1/2, recall that q + = 
1/2. Eq. © then yields 

cosk x + As)n kBo - + XAn^Aa] 

-JvY^ + 3/*+AaV (21) 

Therefore —dF/d/j, = 3N, as is required to ensure three 
electrons per unit cell. Further, one finds that n^Ba = 
for all k and a from the requirement that F is an ex- 
tremum with respect to As , which we expect from qs = 1 
and the constraint of one particle per cite. Finally, 
Ska- n kAa = N due to the requirement that F is an 
extremum with respect to Xa, which is simply the ex- 
pected constraint that the A sublattice is half filled with 
quasiholes. Note, however, that because of the infinite 
U or, equivalently, the constraint of one particle per site, 
the half filled A sublattice is insulating. This is also re- 
flected in the fact that qAA = Qab = 0, i.e. there is zero 
quasiparticle weight for fermions on the A sublattice. So 
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the state qA — 0, qs = 1 (or, equivalently, q + = 1/2, 
g_ = —1/2) describes a charge transfer insulator, which 
has the characteristics of both Mott (A sublattice) and 
band (B sublattice) insulators. Thus, we find that, 



F A 
N = ~2 ~ 



(22) 



which is clearly just the classical energy of paramagnetic 
insulating state. 



Comparison with the empirical theory of Choy 

et al. 




t>0, (n a -r^)/2, exact diag., U=100 
t>0, q_=(n A -n^/2, slave boson, L=100 
t>0, q_=(n A -n^/2, slave boson, L=1000 
t>0, (n -n + )/2, slave boson, L=1000 
t<0, (n -ry/2, exact diag., U=100 
t<0, q_=(n A -n B )/2, slave boson, L=100 
t<0, q =(n A -ry/2, slave boson, L=1000 
t<0, (n -nJ/2, slave boson, L=1000 
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As discussed in section ITBl Choy et al^ have studied 
an empirical model Hamiltonian for holes and doublons 
motivated by the presumed low energy processes about 
an antiferromagnetically ordered charge transfer insulat- 
ing state. This empirical model gives good agreement 
with the measured temperature dependence of both the 
resistivity and the Hall coefficient. It is therefore interest- 
ing to compare our Eq. ^ with the Hamiltonian studied 
by Choy et al. [their Eq. (2)]. 

Upon making the following notational changes to and 
substitutions into Eq. tq A A — > — tfj 4 , tq BB — > 

A- — * V + J; tq AB — > t: fiBa — > h la ; J2a fiAafjAa = 

S ij - KaPjA'i biA -> (k; b' iA fjBo = fiAafjBa and neglect- 
ing constant terms one finds that the two Hamiltonians 
are identical. Notice, in particular, that the renormalised 
ionic potential, A_, in our theory plays the role of the 
inter-site interactions, V + J, in that of Choy et al. Also 
their non-magnetic regime (a = 0) corresponds to q- = 
in our slave boson theory. The ground state studied by 
Choy et al. corresponds to the charge transfer insulator 
we only find for \A/t\ ^> 1. Thus, at least in this regime, 
the two theories appear to be equivalent. In order to 
get agreement with the experimentally measured tem- 
perature dependence of the resistivity and the Hall coef- 
ficient Choy et al. require a specific temperature depen- 
dence of (in our notation) q A A, <7bb, and A_. Although 
we will not present finite temperature calculations below, 
such a temperature dependence is a quite reasonable ex- 
pectation for our mean field theory. Indeed slave boson 
mean-field theories of other strongly correlated models 
have boson mean- fields, similar to g_(T), that decrease 
monotonically with increasing temperature and vanish at 
some coherence temperature, T* -C T F , the Fermi tem- 
perature of the non-interacting system i 54 ' 55 

III. NUMERICAL RESULTS 

We have solved the self consistent equations (fTB|) on an 
L x L reciprocal space mesh for L = 100 and 1000. For 
the L = 1000 solutions we also tightened the convergence 
criteria by a factor of 50. We find that the solutions are 
well converged as the difference between the L = 100 and 



FIG. 2: [Color online] Charge disproportionation as a func- 
tion of |A/t| for both signs of t. Here we compare the real- 
space charge disproportionation, ua~ tib, calculated by exact 
diagonalisation2£ with the results from the slave boson calcu- 
lations where q~ = (ua ~ n_g)/2. The two methods are in 
excellent quantitative agreement for small |A/£|, but at large 
|A/t| the slave boson theory predicts complete charge dispro- 
portionation, i.e., one sublattice becomes completely filled. In 
the exact diagonalisation results complete charge dispropor- 
tionation is only approached asymptotically as |A/tj — > oo. 
In the slave boson theory the system is metallic except when 
there is complete charge disproportionation. In the exact di- 
agonalisation calculations a gap opens when there is only a 
small charge transfer between the two sublattices. We also 
plot the reciprocal space (band) charge disproportionation, 
n_ — n+. We see that the MIT does not occur at the same 
|A/t| as the bonding band becomes filled, in spite of this im- 
plying that the antibonding band is half filled and the fact 
that U — oo. This emphasises the fact that the half filled 
band is not equivalent to a half filled ID chain. We plot q_ 
(similar results are found for other parameters) for two differ- 
ent sized (L x L) reciprocal space meshes. In the L = 1000 
case the convergence criteria is also tightened by a factor of 50, 
which demonstrates that the calculations are well converged. 

L = 1000 self consistent solutions is no larger than a few 
parts in a thousand for any of the parameters. By way 
of an example, in Fig. [2] we report both sets of results, 
it can be seen the difference in the results is significantly 
smaller than the width of the lines. Therefore, we report 
only the L — 1000 results in the subsequent figures. 

A. Charge disproportionation and the 
metal-charge ordered insulator transition 

In Fig. [2] we plot the real space charge disproportion- 
ation, g_ = [n A - tib)/2, where n M = {J2ia 
is the average number of electrons per sublattice site, 
against |A/t|. It can be seen that the results are strongly 
dependent on the sign of t. The system is metallic 
for \q_\ < 1/2; \q_\ = 1/2 corresponds to exactly one 
fermion on each site of the A-sublattice and exactly 
one boson on each site of the B-sublattice. The metal- 
insulator transition (MIT) occurs at \A/t\ — 8.0 for t < 
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and |A/f| = 5.0 for t > 0. 

We also plot the reciprocal space charge disproportion- 
ation (band filling) in Fig. [2l For t > the bonding band 
is filled for all non-zero A/t (for A = the unit cell is 
halved so there is only a single band) . For t < the bond- 
ing band is only partially filled for | A/t \ < 1.5. As we will 
discuss in detail below, the quasiparticles in the bonding 
(-) band are significantly lighter than those in the, nearly 
half-filled, antibonding (+) band. It is interesting to ob- 
serve therefore that for t < 0, 1.5 < \A/t\ < 8.0 and 
t > 0, < \A/t\ < 5.0 the antibonding band is half-filled 
and, even though U = oo, the system remains metallic. 
At first sight this is rather surprising as the system is 
two-dimensional and even quasi-one-dimensional. How- 
ever, because the bands arise from the hybridisation of 
the A and the B sublattices the usual real space argu- 
ments that explain the Mott insulating state do not ap- 
ply here. There is no interaction between the quasipar- 
ticles other than the real space constraint of one fermion 
(or boson) per site. Therefore, unless there is complete 
charge disproportionation the system remains metallic. 

It is interesting to contrast the results of the slave bo- 
son calculations with our previous exact diagonalisation 
results for finite clusters with U — 100|t|i 23 ' 24 At least at 
low | A/t |, the slave boson theory gives excellent agree- 
ment with the exact diagonalisation for some physical 
quantities, e.g., the level of charge disproportionation (cf. 
Fig. [2j. But, there is an important qualitative differ- 
ence between the slave boson calculations and the ex- 
act diagonalisation results: the nature of the insulating 
state. In the slave boson calculations an insulating state 
is only realised when there is complete charge dispropor- 
tionation leaving one sublattice completely filled and the 
other half-filled. The filled sublattice acts as a band in- 
sulator and, because U — oo, the half filled sublattice 
becomes a Mott insulator. This charge transfer insulator 
is very different from the covalent insulating state pre- 
dicted from the exact diagonalisation calculations^^ 
These calculations predict that the insulating state oc- 
curs at quite small values of |A/t| in spite of there be- 
ing only rather weak charge disproportionation. This 
insulating state depends crucially on the hybridisation 
between the two sublattices and is analogous to a cova- 
lent insulator i££ Our mean field theory neglects non-local 
correlations that are included in the exact diagonalisation 
calculations and cause the covalent insulating state. Pre- 
liminary results suggest that DMFT, which also neglects 
non-local correlations but is more sophisticated than our 
mean field theory, cannot describe the covalent insulator 
phase either. 

In Figs. [3] and [4] we plot the ground state energy, F, 
as a function of g_ for various values of |A/t| and the 
self consistently calculated curve of F against with 
|A/t| as a parametric parameter. It can be seen from 
these plots that the self consistent solutions are the true 
minima. Further, the absence of double well structures in 
the plots of F against q- strongly suggest that the MIT is 
a second order phase transition in the slave boson theory. 
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FIG. 3: [Color online] The ground state energy, F, relative to 
that in the insulating state for the same value of A as a func- 
tion of g_ for various values of A and t < 0. The solid black 
line indicates the free energy found at the self consistent solu- 
tion, with crosses marking the data for values of A at which 
we plot the full F(q-) curve. This plot indicates that the self 
consistent solutions are indeed minima of F with respect to 
the variation of q~ [note that Eq. (|16ajl only requires that the 
self consistent solution is a turning point in F with respect to 
q~, whereas, physically, we seek the minimum]. These results 
also strongly suggest that the metal insulator transition in the 
slave boson theory is second order as no double well structure 
is seen in F(q-). 
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FIG. 4: [Color online] Same as Fig. [3]but for t > 0. Again this 
shows that the self consistent solutions are true minima, and 
that the metal-charge ordered insulator transition is second 
order. 



B. Band structure of the fermionic quasiholes 

We now turn to analyse the properties of the fermionic 
quasiholes. As these are non-interacting particles one can 
straightforwardly calculate the properties of the quasi- 
holes. It is interesting to compare these results with the 
band structure of the non-interacting (U = 0) ionic Hub- 
bard model, which we discuss in Appendix [Al 

Before discussing the band structure in detail we 
should point out that the bands are, to some extent, an 
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FIG. 5: [Color online] Band structure of the fermions for 
t < 0. The band structure is not unlike that of the non- 
interacting system (Fig. I14p . However, for the fermions, in- 
creasing \A/t\ has a much stronger effect than it does for the 
non-interacting electrons. In particular as \A/t\ is increased 
the antibonding band (solid lines) is flattened and eventually 
becomes flat when the metal-insulator transition is reached at 
\A/t\ = 8.0. The high symmetry points of the first Brillouin 
□ 
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FIG. 6: [Color online] Band structure of the fermions for t > 
0. There are significant differences from the band structure 
of t < (Fig. [5]l . In particular the metal-insulator transition 
occurs at only |A/t| = 5.0. Also see the online moviei^ 



non-interacting system. It can be seen that at A = 0, 
for both signs of t, the quasiparticle bands are signifi- 
cantly narrower than the equivalent non-interacting elec- 
tron bands. As |A/t| increases the antibonding band 
narrows and eventually becomes flat at the MIT. While 
the bonding band broadens somewhat as |A/t| increases. 

These conclusions are consistent with changes observed 
in the density of states, which we plot in Figs. [7] and [5] 
(also see the animations online^), calculated by treating 
the fermions as non-interacting particles. These figures 
should be contrasted with the equivalent plots of the non- 
interacting system (Figs. IT51 and [TO]) . For both signs of t 
as |A/i| is increased the antibonding band narrows into 
the quasiparticle peak, familiar from the Mott-Hubbard 
transition^ Note that, as the fermions are hole-like, the 
states with E > are filled and the states with E < 
are unoccupied. 

We plot the the Fermi surfaces of the quasiholes in 
Figs. [5] and [TO] Both signs of t exhibit strong nest- 
ing, suggesting that there may be magnetic instabilities, 
which we do not consider here. For t < and |A/t| < 1.5 
there are small two-dimensional Fermi pockets of light 
quasiholes arising from the bonding band, as well as the 
quasi-one-dimensional sheets arising form the antibond- 
ing band. For \A/t\ — > the small Fermi pocket occupies 
0.7% of the Brillouin zone. As \A/t\ increases the Fermi 
pocket shrinks and it vanishes at \A/t\ = 1.5. 

In the A0.5C0O2 materials the unit cells are twice the 
size of that of our model. This results in a zone folding 
of the Brillouin zone in the ^-direction, cf. Fig. [TJ This 
will result in the formation of additional small pockets 
around the X and M points (of the A0.5C0O2 Brillouin 
zone, Fig. [If), which are formed from fermions in the 
antibonding band. If one simply folds our results into the 
A0.5C0O2 Brillouin zone the pockets from the bonding 
band occupy 1.4% of the Brillouin zone for \A/t\ — > (as 
the Brillouin zone has halved in size, while the pocket 
area remains unchanged) and the pockets around the X 
and M points arising from the antibonding band occupy 
8.6% and 9.9% of the Brillouin zone respectively. 



artefact of the mean field approximation. In the exact 
solution correlations that are not captured by the mean 
field approximation may destroy the description of the 
system in terms of extended bands (cf. Ref. I24T ). Nev- 
ertheless, if the mean field theory does provide a rea- 
sonable description of the ionic Hubbard model then one 
expects that some of the features of the band structure 
described below do survive in the exact solution. Fur- 
ther, the bands are well defined in the mean-field slave- 
boson theory. Therefore, it is legitimate to calculate their 
properties in this context and this is extremely useful for 
building intuition about this problem. 

Figs. [5] and [6] show how the dispersion relations of the 
fermionic quasiholes vary with |A/t| for each sign of t 
(also see the animations online^). These should be com- 
pared with the equivalent plots (Figs. [T4"l and [TO]) for the 



C. Spectral density 

The combination of small \ A/t\ and large U/\t\ signifi- 
cantly changes the dispersion relations of the quasiparti- 
cles from those of the non-interacting electrons, compare 
section IIIIBI with Appendix [5] and particularly Figs. [5] 
and [6] with Figs. [14] and [15] First, the antibonding 
band has a much weaker dispersion than for [7 = than 
for U = 00. This is seen in ARPES experiments on 
Nao.5Co02 and K0.5C0O2, where the measured Fermi ve- 
locity is several times smaller than the value calculated 
from the LDA with no sodium ion ordering (cf. Fig. 2i 
of Ref. @). As A increases weight gets shifted from the 
antibonding to the bonding band. 

We plot two different reciprocal space cuts of the 
calculated spectral function for representative values of 
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FIG. 7: [Color online] Density of states (in units where \t\ = a = 1) of the fermions for t < 0. Dashed (blue and pink) lines 
indicate the contributions of the individual bands. The narrowing of the antibonding band, as \A/t\ increases, is clearly visible 
from these plots, while it can be seen that the bonding band actually becomes a little wider as \A/t\ increases. Clearly, the 
DOS is very different from the non-interacting case, cf. Fig. 1181 and these difference become greater as \A/t\ increases. Also 
see the online movie.— 



|A/t| and t < in Figs. QT] and [T2] (also see the ani- 
mations online—). As shown in Fig. [1] the unit cells of 
A0.5C0O2 are twice as large as the unit cell of Hamilto- 
nian In order to facilitate comparison with experi- 
ment, we have plotted these results in the Brillouin zone 
of Nao.5CoC>2 by 'folding over' the results for Hamilto- 
nian j2|) into the Brillouin zone of the actual materials. 
At T = Eq. (|18[) gives a plane of 5-functions (in the 
three dimensional k, u> space). It is well known that 



S(e — uj) = lim 



exp [-(( 



2 /46] 



2V7T& 



(23) 



Therefore we have broadened the results by replacing the 
(^-function with a Gaussian with b = O.Oli 2 . We report 
our results with the Fermi factor appropriate for ksT = 
|*|/60 ~ 20 K for -t = 0.1 eV. 

It can be seen from Figs. ITTI and [T2l that as \A/t\ is 
increased the antibonding band remains at the Fermi en- 
ergy, but loses spectral weight, consistent with the cal- 
culated density of states for fermions (Figs. [7] and EJ. 
Meanwhile, the bonding band is shifted down in energy. 
In the T-M cut (Fig. [TTj) we see that the antibonding 
band narrows and loses intensity as \A/t\ is increased. 



In the T-K cut (Fig. IT2")) we see that the bonding band 
moves from the Fermi level at A = to lower energies 
such that it is completely filled at larger values of |A/i|. 
At first sight Fig. [T^] might appear to suggest that the 
bonding band also narrows as |A/£| increases, however it 
is clear from Fig. QT] that this is not the case. Rather, 
as \A/t\ increases the system becomes more one dimen- 
sional (for t/A = 0, the system consists of uncoupled 
one-dimensional chains), increasing the range of in 
the cuts shown in Fig. [11] and decreasing the range of 
in the cuts shown in Fig. [12] (cf. Figs. [5] and [5]). 



D. Quantum oscillations 

As we are considering a quasi-two-dimensional metal 
the cyclotron effective mass, m±, associated with 
Shubnikov-de Haas oscillations of the part of the Fermi 
surface arising from the ± band is^ 



m± = 2Trh 2 D±(E F ), 



(24) 



where D±(Ep) is the renormalised density of states per 
spin at the Fermi energy, Ep, arising from the ± band. 
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FIG. 8: [Color online] Density of states (in units where \t\ — a = 1) of the fermions for t > 0. Dashed (blue and pink) 
lines indicate the contributions of the individual bands. As in Fig. [7] we see that the antibonding band rapidly narrows at 
\A/t\ increases. Furthermore, the width of the bonding band is almost independent of |A/i|. Again, the DOS is very different 
from the non-interacting case, cf. Fig. 1191 particularly for large |A/t|. Also see the online movie.— 
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FIG. 9: [Color online] Fermi surface of the fermions for t < 0. 
For small |A/t| there are small Fermi pockets around the Y 
points. Comparing the Fermi surface with the band struc- 
ture of the fermions shows that these pockets arise from the 
bonding band. One therefore expects that these fermions are 
lighter than those in the quasi-one-dimensional sheets, which 
arise from the antibonding band (cf. Fig. [13) . The Brillouin 
zone is that of the ionic Hubbard model (cf. Fig. fjj). 



We plot the variation of the cyclotron mass with \A/t\ for 
parameters relevant to Nao.sCoC^, i.e., t = —0.1 eV, a 
basal lattice constant, a = 2.82 A, and a unit cell basal 
area of \/3a x 2a, in Fig. [13] The cyclotron mass for the 
bonding band associated with the small pockets is about 
one quarter of that for heavy electrons associated with 
the antibonding band. 

In order to observe the states arising from the anti- 
bonding band, i.e., the quasi-lD sheets on the Fermi sur- 
face, in quantum oscillation experiments a field that is 
sufficiently large to cause magnetic breakdown between 
the two bands must be applied. The field required for 
magnetic breakdown and the large effective mass associ- 
ated the antibonding band make it more difficult to ob- 
serve the breakdown orbit than than the small pockets 
in quantum oscillation experiments. 

For realistic parameters for Nao.sCoC^, A ~ \t\ <C U 
and t < (cf . section II Bl) , the slave boson theory pre- 
sented above suggests that the system is metallic. How- 
ever, an important question, which we cannot conclu- 
sively answer because of the uncertainty in the precise 
value of |A/f| relevant to Nao.sCoC^, is: is |A/i| large 
enough in Nao.sCoC^ that the bonding band is com- 
pletely filled? For \A/t\ < 1.5 there are two features 
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FIG. 10: [Color online] Fermi surface of the fermions for t > 0. 
For t > the bonding band is always completely filled. There- 
fore this Fermi surface arises entirely from the antibonding 
band. As \A/t\ is increased only rather subtle changes are 
observed in the Fermi surface, most notably, the nesting is 
slightly decreased. The Brillouin zone is that of the ionic 
Hubbard model (cf. Fig. []J) . 

in the Fermi surface (cf. Fig. [9]): small hole- like pock- 
ets of light fermions arising from the bonding (nearly 
filled) band, and quasi-one-dimensional sheets of heavier 
electron-like quasiparticles. This seems to be quite con- 
sistent with quantum oscillations, 7 which are observed 
to have a frequency correspond to orbits that encompass 
less than 1% of the first Brillouin zone, in apparent vio- 
lation of Luttinger's theorem. A natural explanation of 
this experiment is that Balicas et alX are only able to ob- 
serve the, small, light, quasiparticles in the bonding band 
and do not see oscillations arising from the antibonding 
band, either because the quasiparticles in the antibond- 
ing band are two heavy, or because the field they applied 
is insufficient to cause magnetic breakdown of the 'gap' 
between the quasi-one-dimensional sheets and the two- 
dimensional pockets. However, for \A/t\ > 1.5 the small 
Fermi pockets are absent as the bonding band is filled, 
and our results would be inconsistent with the quantum 
oscillations observed by Balicas et al. 



E. A picture of Nao.sCoCh 

The experimental situation, particularly in Nao.5Co02, 
is somewhat confusing. Several experiments suggest an 
insulating state: resistivity versus temperature has a 
negative gradient at low temperatures^ and angle re- 
solved photoemission spectroscopy (ARPES) 6 and opti- 
cal conductivity 5 - both observe a gap, all of which sug- 
gests that below 51 K Nao.sCoC^ is a insulator. However, 
quantum oscillations are observed) 7 - which suggest that 
there are quasiparticles and that the system is metallic 
with a small Fermi surface (occupying less that 1% of 
the cross-sectional area of the first Brillouin zone) and 
an effective mass similar to the free electron mass. 



Below we propose a picture, inspired by both the mean 
field slave boson theory (above) and our previous exact 
diagonalisation calculations for the same ionic Hubbard 
model ) 23 ' 24 that offers a possible reconciliation of these 
seemly contradictory experiments. However, it is impor- 
tant to stress that this is a tentative proposal and does 
not correspond to the solution found in either theory. 
Nevertheless, we will argue, below, that this picture is 
plausible given the limitations of the mean field slave 
boson theory and the finite size effects in the exact di- 
agonalization calculations. Further calculations will be 
required to discover whether this proposal is indeed sup- 
ported by the theory of the ionic Hubbard model. 

The fundamental definition of a metal is a substance 
with charge carrying excitations at arbitrarily low ener- 
gies above the ground state. In most situations it follows 
that the resistivity monotonically increases with temper- 
ature as increasing the temperature decreases the quasi- 
particle lifetime. However, this is not necessarily the 
case. For example, in many strongly correlated systems, 
such as the organic superconductors^ the charge trans- 
port is incoherent at high temperatures - this leads to a 
resistivity with a broad maximum slightly above the tem- 
perature where quasiparticle coherence is losti 6 - 1 - Above, 
we have shown that the slave boson theory of Hamil- 
tonian Q produces two bands. We will assume some 
features of these bands survive even in the exact solu- 
tion even if the simple single particle states, of the slave 
boson theory, do not. The antibonding band is narrow 
and it is likely that coherent quasiparticles do not emerge 
from these states until quite low temperatures. Within 
slave boson theories the coherence temperature can be 
estimated from the Bose condensation temperature, Tb 
of the slave bosons^ 2 - which, in 2D, is proportional to 
the density of bosons, ps- For |A/t| — » there are 0.993 
quasiholes per unit cell in the antibonding band, corre- 
sponding to pb = 0.007, and the density of bosons de- 
creases as \A/t\ is increased, going to zero at |A/i| = 1.5. 
Estimating the Bose condensation temperature (cf. Refs. 
[52l and l62l) by either Tb ~ (27rps)/(v / 3a 2 ms), where a is 
the basal plane lattice constant (2.82 A for Nao.sCoC^) 
and tub ~ 2m is the mass of the boson (which is the sim- 
plest approximation given that the boson corresponds to 
two electrons), or Tb ~ iirpst yields Tb ~ 100 K for 
Pb = 0.007. However it is well known that these simple 
formulae dramatically overestimate the coherence tem- 
perature as they neglect gauge fluctuations . 52 ' 62 Thus the 
coherence temperature (T* <C 100 K) is expected to be 
extremely low in the experimentally relevant parameter 
regime. 

This low coherence temperature is consistent with the 
observed incoherent transport evident from the resistiv- 
ity above ~51 K. As the temperature is lowered below 
51 K the resistivity rises rapidly, suggesting that a gap 
opens. The slave boson calculations do not predict a 
true gap unless |A/t| 3> 1. Experiments observe mag- 
netic ordering , 26 i 31 i 37 i 41 which we have not considered 
here, below 88 K, but, as the gap does not appear to 
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A = 0;i<0 A = 2jtj; t < A = 4|t|; t < 




FIG. 11: [Color online] Cut of the spectral density, A(\t,uj), along the F-M direction of the triangular lattice Brillouin zone 
(cf. Fig. QJl). A(\i,u)) is calculated from Eq. (|18|l . broadened as described in the text below Eq. (|23[1 . and plotted for 
ksT — |t|/60 ~ 20 K. Green dotted lines are the slave boson dispersion, blue dashed lines are the bare dispersion. The abscissa 
is in A - , with the lattice constant taken as a = 2.82 A, as is appropriate for Nao.sCo02.— As |A/t| is increased the spectral 
weight is transferred from the antibonding band to the bonding band and the antibonding band narrows. Also see the online 
movie.— 



open until 51 K, one cannot conclude that the magnetic 
order causes the gap to open. However, previous exact 
diagonalisation calculation o 23 ' 24 suggest that in the ex- 
perimentally relevant parameter regime, — t ~ A <C U, 
non-local correlations, not captured by our mean field 
treatment, drive the formation of a covalent insulator. 
We propose that a gap opens on the antibonding band 
at 51 K due to non-local correlations, analogous to the 
covalent insulator. But, we also propose that the bonding 
band remains ungapped. 

This proposal is rather similar to the theory of 
the orbital-selective Mott transition. The question of 
whether, in multi-band systems, one band can becoming 
insulating, due to electronic correlations, while others re- 
main metallic has been extensively discussed recently, 63 
motivated, in part, by the volume collapse in Ce (Ref. 

the heavy fermion behaviour of Ca 2 _ x Sr x Ru04 
(Ref. f65h and the nodal-anti-nodal differentiation in the 
cuprates^ 6 - Although the debate is ongoing, the current 
theoretical consensus is that orbital-selective Mott tran- 
sitions can occur for some parameter regimes of some 
models.— Therefore, testing our proposed picture will 
require calculations beyond the mean-field slave-boson 
theory presented above. 

Given the important differences between the current 
mean field calculations and our previous exact diago- 
nalisation calculations, particularly the critical value of 
|A/t| for the metal-insulator transition, it is important 



to ask, whether this explanation is likely to hold in the 
exact solution and, more importantly, for the real mate- 
rial. Firstly, we note that the real space charge dispro- 
portionation calculated from our mean field theory is in 
very good agreement with exact diagonalisation calcula- 
tions in the experimentally relevant regime. This gives 
one hope that the filling of the bands calculated from 
the mean field theory may also be reasonable (although 
non-local Coulomb correlations could cause additional 
hybridisation effects, cf. Ref. [24]). Secondly, we note 
that the finite size clusters that were, necessarily, studied 
in the exact diagonalisation calculations mean that it is 
extremely difficult to accurately study the bonding band. 
Taking a, typical, 18 site (9 unit cell) cluster, a density of 
0.007 holes per unit cell means that one would expect less 
than 1 / 16 th of a hole in the bonding band in the exact 
diagonalisation calculation. However, exact diagonalisa- 
tion calculations do accurately reproduce the band filling 
of the analytical exact solution for U — 0. 23 <i2i Therefore, 
it is not clear, at present, whether the exact diagonali- 
sation calculations can resolve the tiny hole pockets sug- 
gested by our calculations and the quantum oscillation 
experiments. 

As the exact diagonalisation calculations do not rule 
out the possibility that small weakly correlated Fermi 
pockets containing a tiny number of charge carries re- 
main even in the low temperature 'insulating' state, let 
us briefly discuss how this would produce a consistent 
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A = 0;i<0 A = 2jtj; t < A = 4|t|; t < 




FIG. 12: [Color online] Cut of the spectral density, A(k,u>), along the T-K direction of the triangular lattice Brillouin zone 
(cf. Fig. QJl). A(\i,u)) is calculated from Eq. (JTSJ, broadened as described in the text below Eq. (|23[1 . and plotted for 
ksT — |t|/60 ~ 20 K. Green dotted lines are the slave boson dispersion, blue dashed lines are the bare dispersion. The abscissa 
is in A - , with the lattice constant taken as a = 2.82 A, as is appropriate for Nao.sCo02.— At small [A/t| the bonding band 
crosses the Fermi energy in this direction (cf. Figs. [5] and [9}. As \A/t\ is increased the bonding band is pushed below the 
Fermi energy and spectral weight is transferred from the antibonding band to the bonding band, which is not dispersive in this 
direction. Also see the online movie.— 



explanation of the experiments on Nao.sCoC^. 

It is important to note that the two bands act as par- 
allel channels for charge transport. Above 51 K charge 
transport in the antibonding band will be incoherent (as- 
suming the coherence temperature < 50 K, as we have 
argued above) and the resistivity in this channel can be 
estimated to be > hc/e 2 , where c is the interlayer spac- 
ing, i.e., at or above the Mott-Ioffc-Rcgcl limit. Below 51 
K a gap opens on the antibonding band, but as the gap is 
small, 0(t) (Ref. [23,24) the resistivity in the antibonding 
channel is not increased by orders of magnitude from the 
resistivity value in the incoherent metallic state. How- 
ever, because of the extremely low charge carrier concen- 
tration in the bonding band (metallic channel) the con- 
tribution of the bonding band to the conductivity will be 
even smaller than that of the antibonding band (incoher- 
ent/insulating channel). Therefore the metallic channel 
will be unable to 'short circuit' the incoherent/insulating 
channel. This contention is supported by the fact that 
even above 51 K the resistivity is larger than the Mott- 
Ioffe-Rcgel limit. Thus the observation that the resistiv- 
ity decreases with increasing temperature is not neces- 
sarily inconsistent with the observation of quasiparticles 
in quantum oscillation experiments. 

Foo et alA found that the conductivity, a, of 
Nao.68CoC>2 is very similar to that of Nao.3iCo02 and 
that both of these are an order of magnitude larger than 



the conductivity of Nao.sCoC^. As a cx n, one expects 
that the contribution of the electrons in the bonding 
band to the conductivity of Nao.sCoC^ will be two or- 
ders of magnitude smaller than the total conductivity of 
Nao.68Co02 or Nao.3iCo02. Therefore Foo et al.'s results 
are consistent with the, counterintuitive, picture that the 
insulating band short-circuits the metallic band! 

ARPES, which measures A(k, u>), has been reported in 
the directions T-M and T-K, i.e., in the directions of the 
cuts shown in Figs. 1111 and [T2l respectively. However, our 
slave boson calculations do not capture the gap on the 
antibonding band required in our proposal. Therefore the 
gap observed in the ARPES is likely to correspond to the 
'covalent insulator' gap on the antibonding band. Opti- 
cal conductivity does not observe a Drude peak, however 
the height of the Drude peak is proportional to the charge 
carrier density, which would make the tiny Fermi pock- 
ets suggested by our slave boson calculations difficult to 
observe. 

As one expects that the Fermi pockets do not con- 
tribute significantly to cither the resistivity or the Hall 
coefficient the correspondence of our mean-field slave bo- 
son theory to the empirical low energy Hamiltonian stud- 
ied by Choy et alM (discussed in section lHCj) become im- 
portant. The non-local correlations on the antibonding 
band open up a gap on the quasi-one-dimensional sheets 
of the Fermi surface. This provides the gap assumed by 
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FIG. 13: [Color online] Variation of the renormalised cy- 
clotron effective mass, m±, in units of the free electron mass, 
m, with |A/t| for t < 0. Note that in the main panel the ordi- 
nate is plotted on a logarithmic scale: the same data is plotted 
on a linear scale in the inset. m_/m ~ 1 and m+/m ~ 4 for 
small |A/t|. The cyclotron effective mass of the antibonding 
band diverges as the metal-charge transfer insulator transi- 
tion is approached, consistent with our finding that the metal- 
charge ordered insulator transition is second order, cf. Fig. [3] 
The cyclotron effective masses are calculated from Eq. (|24p 
with parameters appropriate for Nao.sCoO^: t = —0.1 eV, an 
in-plane lattice constant, a — 2.82 A, and a unit cell basal 
area of y/3a x 2a. 



Choy et al. and suggests that their results will carry 
across to our model. This suggests that the measured 
temperature dependence of both the resistivity and the 
Hall coefficient in Nao.sCoC^ have a natural explanation 
in terms of the ionic Hubbard model. However, clearly 
further calculations are require to test the validity of this 
ad hoc proposal and its relevance to Nao.sCoC^. 



FIG. 14: Band structure of the ionic Hubbard model for 
U = 0, at three-quarters filling, and t < 0. Increasing 
|A/£| makes the system increasingly anisotropic, which serves 
to flatten the bands slightly in some directions. The con- 
straint of three quarters filling 'pins' the antibonding band 
to the Fermi energy, but increasing [A/i| pushes the bonding 
band to lower energies. 



low density of charge carriers in the bonding band mean 
that the metallic bonding band is not able to 'short cir- 
cuit' the incoherent/insulating antibonding band. Cal- 
culations beyond the mean field slave boson calculations, 
which we have presented above, will be required to test 
this picture. 

Another important avenue for future work will be to 
allow for magnetic ordering in the theory as this may 
improve the degree of agreement with experiments and 
allow one to compare with other experiments such as 
NMRM and neutron scattering^! 



IV. SUMMARY AND CONCLUSIONS 



We have presented a mean-field slave boson theory of 
the ionic Hubbard model on the triangular lattice with 
alternating stripes of site energy. This model has two 
bands: one of which is weakly correlated and nearly filled 
or filled, the other is nearly half-filled or half-filled and 
hence strongly correlated. The results depend strongly 
on the sign of t. The light band is always filled for t > 0, 
but only becomes filled at |A/f| = 1.5 for t < 0. A 
metal-charge transfer insulator transition occurs at larger 
|A/i| (5.0 for t > and 8.0 for t < 0), when complete 
charge disproportionation occurs and one sublattice is 
filled and the other in half filled. 

We have also proposed a speculative picture of 
Nao.5Co02- In particular, we have argued that the ob- 
served quantum oscillations arise from quasiparticles in 
the bonding band, but these are not seen in the zero 
field resistivity and several other experiments because the 
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FIG. 15: Band structure of the ionic Hubbard model for U = 
0, at three-quarters filling, and t > 0. This band structure 
is markedly different from the t < case (Fig. 1 14f) . This 
is consistent with the general expectation that particle-hole 
symmetry is absent on frustrated lattices and hence that the 
sign of t has strong effects on the physics. 



APPENDIX A: NON-INTERACTING 
ELECTRONS 

For [7 = there are two bands, bonding (-) and anti- 
bonding (+), and whose dispersion relations are 



-2t cos(k x ) 



±^A 2 /A+(4tcos(k x /2)cos(k y /2)) 2 , (Al) 

where k x and k y are defined in the reduced (1 x -\/3) Bril- 
louin zone of Hamiltonian @, shown in Fig. [T^. These 
dispersion relations are plotted in Figs. [T3] and [15] for 
t < and t > respectively. 

To calculate where the bonding band becomes filled we 
need to consider the maximum of £^ . This is complicated 
by the fact the chemical potential, /i, is a function of A. 
For t > the top of the bonding band is at = = 
eV s. Thus we find that the bonding band become full 
at A > A> = 2^(A > )+4|*|. However, for A = 0+ or (T 
and t > one finds that /.i = — 2.0t and this condition is 
satisfied. Therefore the bonding band is filled, and hence 
the antibonding band is half filled, for an infinitesimal 
A. Note that this is not the case for A = as we require 
at least an infinitesimal A to double the size of the unit 
cell and cause the splitting into two bands, cf. Figs. [T]a, 
and [T|d. For t < the top of the bonding band is £y 
and we find the bonding band is filled for A > A < = 
2/i(A < ) + 4|i|. Solving this condition numerically we find 
that A< = -0.64|t|. 

In Figs. [TBI and [T7l we plot the Fermi surfaces for non- 
interacting electrons for each sign of t at a range of |A/t|. 
The first point to note is that at low |A/i| the two Fermi 
surfaces are very different. As |A/i| — * oo both Fermi 



FIG. 16: [Color online] Fermi surface for non-interacting elec- 
trons for t < 0. For small \A/t\ there are small Fermi pockets 
around the Y points. As \A/t\ is increased the Fermi surface 
becomes increasingly quasi-one-dimensional. 
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FIG. 17: [Color online] Fermi surface for non-interacting elec- 
trons for t > 0. For t > the bonding band is always com- 
pletely filled. Therefore this Fermi surface arises entirely from 
the antibonding band. Note that the Fermi surface is very dif- 
ferent from that for t < (Fig. [15) . 



surfaces become quasi-one-dimensional; but, even in this 
limit the two Fermi surfaces are different in important 
ways as for t < the Fermi surface is electron like, while 
of t > the Fermi surface is hole like. Only when we 
reach the limit, and both Fermi surfaces are straight lines, 
are the two Fermi surfaces the same. 

We plot the DOS of the non-interacting (U — 0) system 
in Figs. [18] and [19] for t < and t > respectively. Notice 
that the plots are very different for the two different signs 
of t. In particular at A = the Fermi energy is at the 
van Hove singularity for t > 0, but is not for t < 0. This 
suggests that the Fermi liquid will be less stable for t > 
than for t < 0. Indeed we find that, in the slave boson 
theory, the MIT occurs at \A/t\ = 5.0 for t > and 
\A/t\ = 8.0 for t < 0, consistent with this expectation 
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(cf. Fig. [2]). However, as |A/t| is increased the van Hove t > 0, so this, weak coupling, description is clearly not 
singularity is moved away from the Fermi energy even for the whole story of the strong coupling theory. 



* Electronic address: bjpowell@gmail.com 

1 See, for example, N. P. Ong and R. J. Cava, Science 305, 
52 (2004). 

2 J. Merino, B. J. Powell, and R. H. McKenzie, Phys. Rev. 
B 73, 235107 (2006). 

3 K. Takada et al , Nature 422, 53 (2003). 

4 M. L. Foo, Y. Wang, S. Watauchi, H.W. Zandbergen, T. 
He, R. J. Cava, and N. P. Ong, Phys. Rev. Lett. 92, 247001 
(2004). 

5 N.L. Wang et al, Phys. Rev. Lett. 93, 147403 (2004); J. 
Hwang et al, Phys. Rev. B 72, 024549 (2005); S. Lupi et 
al, Phys. Rev. B 72, 024550 (2005). 

6 D. Qian, L. Wray, D. Hsieh, D. Wu, J. L. Luo, N. L. Wang, 
A. Kuprin, A. Fedorov, R. J. Cava, L. Viciu, and M. Z. 
Hasan, Phys. Rev. Lett. 96, 046407 (2006). 

7 L. Balicas, M. Abdel-Jawad, N. E. Hussey, F. C. Chou, 
and P. A. Lee, Phys. Rev. Lett. 94, 236402 (2005). 

8 M. Z. Hasan, Y.-D. Chuang, D. Qian, Y. W. Li, Y. Kong, 
A. P. Kuprin, A. V. Fedorov, R. Kimmerling, E. Roten- 
berg, K. Rossnagel, Z. Hussain, H. Koh, N. S. Rogado, 
M. L. Foo, and R. J. Cava, Phys. Rev. Lett. 92, 246402 
(2004); H.-B. Yang, S.-C. Wang, A. K. P. Sekharan, H. 
Matsui, S. Souma, T. Sato, T. Takahashi, T. Takeuchi, J. 
C. Campuzano, R. Jin, B. C. Sales, D. Mandrus, Z. Wang, 
and H. Ding, ibid. 92, 246403 (2004). 

9 P. Zhang, R. B. Capaz, M. L. Cohen, S. G. Louie, Phys. 
Rev. B 71, 153102 (2005). 

10 M. Roger, D. J. P. Morris, D. A. Tennant, M. J. Gutmann, 
J. P. Goff, J. U. Hoffmann, R. Feyerherm, E. Dudzik, D. 
Prabhakaran, A. T. Boothroyd, N. Shannon, B. Lake, and 
P. P. Deen, Nature 445, 631 (2007). 

11 H. W. Zandbergen, M. L. Foo, Q. Xu, V. Kumar, and R. 
J. Cava, Phys. Rev. B 70, 024101 (2004). 

12 T. Egami, S. Ishihara, and M. Tachiki, Science 261, 1307 
(1993). 

13 A. P. Kampf, M. Sekania, G. I. Japaridze, and Ph. Brune, 
J. Phys.: Condens. Matter 15 5895 (2003); also see ref- 
erences therein and citations thereof for further details of 
the extensive literature on the ionic Hubbard model in one- 
dimension. 

14 A. Garg, H. R. Krishnamurthy, and M. Randeria, Phys. 
Rev. Lett. 97, 046403 (2006). 

15 L. Craco, P. Lombardo, R. Hayn, G. I. Japaridze, and E. 
Miiller-Hartmann, Phys. Rev. B 78, 075121 (2008). 

16 S. S. Kancharla and E. Dagotto, Phys. Rev. Lett. 98, 
016402 (2007). 

17 N. Paris, K. Bouadim, F. Hebert, G. G. Batrouni, and R. 
T. Scalettar, Phys. Rev. Lett. 98, 046403 (2007). 

18 K. Bouadim, N. Paris, F. Hebert, G. G. Batrouni, and R. 
T. Scalettar, Phys. Rev. B 76, 085112 (2007). 

19 K. Penc, H. Shiba, F. Mila, and T. Tsukagoshi, Phys. Rev. 
B 54, 4056 (1996). 

20 G. I. Japaridze, R. Hayn, P. Lombardo, E. Miiller- 
Hartmann, Phys. Rev. B 75, 245122 (2007). 

21 K. Byczuk, M. Ulmke, and D. Vollhardt, Phys. Rev. Lett. 
90, 196403 (2003). 

22 C. A. Marianetti and G. Kotliar, Phys. Rev. Lett. 98 



176405 (2007). 

J. Merino, B. J. Powell, and R. H. McKenzie, Phys. Rev. 
B 79, 161103(R) (2009). 

J. Merino, R. H. McKenzie, and B. J. Powell, 
larXiv:0904.3402l 

F.L. Ning et al, Phys. Rev. Lett. 93, 237201 (2004). 
J. Bobroff et al, Phys. Rev. Lett. 96, 107201 (2006). 
A. J. Williams et al, Phys. Rev. B 73, 134401 (2006); D. 
Argyriou et al, Phys. Rev. B 76, 134506 (2007). 

F. L. Ning et al, Phys. Rev. Lett. 100, 086405 (2008). 

G. Garbarino et al, Rev. B 77, 064105 (2008). 

The dependence on the direction of the field may be ex- 
plained by (i) the larger g-factor for fields parallel to the 
layers, and/or (ii) that the transitions at 51 K and 88 K 
have a larger effect on the bulk magnetic susceptibility par- 
allel to the layers, Yqt) than on Xc (see Fig. 1(c) in Ref. [3ll 
and Fig. 8 in Ref. |32| ) . 

G. Gasparovic et al, Phys. Rev. Lett. 96, 046403 (2006). 

H. Watanabe et al, J. Phys. Soc. Jap. 75, 034716 (2006). 
P. Mendels, D. Bono, J. Bobroff, G. Collin, D. Colson, 
N. Blanchard, H. Alloul, I. Mukhamedshin, F. Bert, A. 
Amato, and A. D. Hillier, Phys. Rev. Lett. 94, 136403 
(2005). 

I. R. Mukhamedshin, H. Alloul , G. Collin, and N. Blan- 
chard, Phys. Rev. Lett. 94, 247602 (2005). 

M.-H. Julien, C. de Vaulx, H. Mayaffre, C. Berthier, M. 
Horvatic, V. Simonet, J. Wooldridge, G. Balakrishnan, M. 
R. Lees, D. P. Chen, C. T.Lin, and P. Lejay, Phys. Rev. 
Lett. 100, 096405 (2008). 

H. Alloul , I.R. Mukhamedshin, G. Collin, and N. Blan- 
chard, Europhys. Lett. 82, 17002 (2008). 
M. Yokoi, T. Moyoshi, Y. Kobayashi, M. Soda, Y. Yasui, 
M. Sato, and K. Kakurai, J. Phys. Soc. Jpn. 74, 1800 

(2005) . 

H. Alloul, I.R. Mukhamedshin, TA. Platova, and A.V. 
Dooglav, Europhys. Lett. 85, 47006 (2009). 
G. Lang, J. Bobroff, H. Alloul, G. Collin, and N. Blan- 
chard, Phys. Rev. B 78, 155116 (2008). 
J. Bobroff et al, Phys. Rev. B. 76 100407(R) (2007). 
T. F. Schulze et al, Phys. Rev. Lett. 100, 026407 (2008). 
T.-P. Choy, D. Galanakis, and P. Phillips, Phys. Rev. B 
75, 073103 (2007). 

K.-W. Lee, J. Kunes, P. Noval, and W. E. Pickett, Phys. 
Rev. Lett. 94, 026403 (2005). 

K.-W. Lee and W. E. Pickett, Phys. Rev. Lett. 96, 096403 

(2006) . 

S. Zhou and Z. Wang, Phys. Rev. Lett. 98, 226402 (2007). 
K.-W. Lee and W. E. Pickett, Phys. Rev. B 76, 134510 

(2007) . 

S. Landron and M.-B. Lepetit. [arXiv:0805.312T1 

D. J. Singh, Phys. Rev. B 61, 13397 (2000). 

C. de Vaulx, M.-H. Julien, C. Berthier, S. Hebert, V. Pra- 

long, and A. Maignan, Phys. Rev. Lett. 98, 246402 (2007). 

M. Capone, L. Capriotti, F. Becca, and S. Caprara, Phys. 

Rev. B 63, 085104 (2001). 

S. E. Barnes, J. Phys. F 6, 1375 (1976); ibid. 7, 2637 
(1977); P. Coleman, Phys. Rev. B 29, 3035 (1984). 



18 




IA/tl=7.9, 1<0, U=0 



2 - 




FIG. 18: [Color online] Density of states for U = and t < 0. Dashed (blue and pink) lines indicate the contributions of the 
individual bands. 



P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 
78, 17 (2006). 

See, for example, section 4.10 of M. L. Boas Mathematical 
methods in the physical sciences (2nd edition, Wiley, New 
York, 1983). 

J. Spalek and W. Wojcik, in Spectroscopy of the Mott 
Insulators and Correlated Metals, Springer Series in Solid 
State Sciences Vol. 119 (Springer- Verlag, Berlin, 1995), pp. 
41-65. 

S. Burdin, A. Georges, and D. R. Grempel, Phys. Rev. 
Lett. 85, 1048 (2000). 

D. D. Sarma, J. Sol. State Chem. 88, 45 (1990). 
Animated versions of Figs. O El El EI 
TT1 and [12] can be viewed, respectively, at 
http : //www . youtube . com/watch?v=818YVFHHs74 , 
http : //www . youtube . com/watch?v=oqbFsEimi3s , 
http : //www . youtube . com/watch?v=q5bmSS3VnaQ , 
http : // www . youtube . com/ wat ch?v=Fz9XKbLZbkE , 
http : //www . youtube . com/watch?v=EjNVEweoOoA , and 
http : //www . youtube . com/watch?v=P2ytkyuzFcA , The 
animations show the same data as is reported in the 
associated figure; each frame of the movie represents a 
different value of |A/t|, specifically each subsequent frame 
represents an increment of 0.1 in |A/t|. 
G. Kotliar and D. Vollhardt, Phys. Today 57, 53 (2004). 
J. Merino and R.H. McKenzie, Phys. Rev. B 62, 2416 
(2000). 

B. J. Powell and R. H. McKenzie, J. Phys.: Condens. Mat- 
ter, 18, R827, (2006). 



J. Merino and R. H. McKenzie, Phys. Rev. B 61, 7996 
(2000). 

P. A. Lee and N. Nagaosa, Phys. Rev. B 46, 5621 (1992). 

A. Liebsch, Europhys. Lett. 63, 97 (2003); Phys. Rev. Lett. 
91, 226401 (2003); Phys. Rev. B 70, 165103 (2004); Phys. 
Rev. Lett. 95, 116402 (2005); A. Koga, N. Kawakami, 
T. M. Rice, and M. Sigrist, Phys. Rev. Lett. 92, 216402 

(2004) ; R. Arita and K. Held, Phys. Rev. B 72, 201102(R) 

(2005) ; M. Ferrero et al., Phys. Rev. B 72, 205126 (2005); 
L. de'Medici, A. Georges, and S. Biermann, Phys. Rev. B 

72, 205124 (2005); M. Ferrero, F. Becca, M. Fabrizio, and 
M. Capone, Phys. Rev. B 72, 205126 (2005); A. Riiegg, 
M. Inderland, S. Pilgram, and M. Sigrist, Eur. Phys. J. 

B. 48, 55 (2005); K. Inaba and A. Koga, Phys. Rev. B 

73, 155106 (2006); P. G. J. van Dongen, C. Knecht, and 
N. Bliimer, Phys. Status Solidi (b) 243, 116 (2006); T. A. 
Costi and A. Liebsch, Phys. Rev. Lett. 99, 236404 (2007); 
N. Bliimer et al., J. Magn. Magn. Mater. 310, 922 (2007); 
K. Bouadim, G. G. Batrouni, and R. T. Scalettar, Phys. 
Rev. Lett. 102, 226402 (2009). 

See, for example, K. Held, A. K. McMahan, and R. T. 
Scalettar, Phys. Rev. Lett. 87, 276404 (2001). 
V. I. Anisimov, I. A. Nekrasov, D. E. Kondakov, T. M. 
Rice, and M. Sigrist, Eur. Phys. J. B. 25, 191 (2002); S. 
Nakatsuji and Y. Maeno Phys. Rev. Lett. 84 2666 (2000). 
M. Ferrero, P. S. Cornaglia, L. De Leo, O. Parcollet, G. 
Kotliar, A. Georges, Phys. Rev. B 80, 064501 (2009). 



19 




E 



FIG. 19: [Color online] Density of states for U = and t > 0. Dashed (blue and pink) lines indicate the contributions of the 
individual bands. 
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